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Abstract. After briefly reviewing recent progress by the HotQCD collaboration in studying 
the 2+1 flavor QCD equation of state, we will focus on results on fluctuations of conserved 
charges by the BNL-Bielefeld and HotQCD collaborations. Higher order cumulants of the net- 
charge distributions are increasingly dominated by a universal scaling behavior, which arises 
due to a critical point of QCD in the chiral limit. Considering cumulants up to the G*** 
order, we observe that they generically behave as expected from universal scaling laws, which 
is quite different from cumulants calculated within the hadron resonance gas model. Taking 
ratios of these cumulants, we obtain volume independent results that can be compared to the 
experimental measurements. We will argue that the freeze-out chemical potentials and the 
freeze-out temperature, usually obtained by a HRG model fit to the measured hadronic yields, 
can also be obtained in a model independent way from ab-initio lattice QCD calculations 
by utilizing observables related to conserved charge fluctuations. Further, we will show that 
the freeze-out strangeness and electric charge chemical potentials can be fixed by imposing 
strangeness neutrality and isospin asymmetry constraints in the lattice QCD calculations, in 
order to accommodate conditions met in heavy ion collisions. All results have been obtained 
with the highly improved staggered quark action (HISQ) and almost physical quark masses on 
lattices with temporal extent of A'^t = 6, 8, 10, 12. 



1. Introduction 

Calculating the phase diagram of strongly interacting matter is one of the most important and 
outstanding problems of non-perturbative QCD. A generic plot of the QCD phase diagram based 
on model calculations and model independent symmetry arguments is shown in Fig{T] To clarify 
whether or not there is a critical point in the phase diagram is relevant for cosmology, astro 
physics and heavy ion collisions. Currently various experimental programs aim at probing the 
phase diagram at nonzero baryon number density, such as the beam energy scan at RHIC, NA61 
at the SPS as well as future experiments at FAIR and NICA. In Fig{T] we have indicated the 
approximate initial conditions of the fireballs as they are created in some of those experiments. 
Also shown is a sketch of the trajectories on which the fireball evolves in the temperature (T) and 
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Figure 1. Generic phase diagram of QCD, based on model calculations and model independent 
symmetry arguments. Also indicated are the approximate initial conditions of the fireballs that 
are created in heavy ion experiments as well as their hydrodynamic trajectories on which they 
evolve, until they reach chemical freeze-out. 

baryon chemical potential (/Us) plane, during its expansion and cooling. It has been shown that 
this evolution can be very well described by hydrodynamics, with a surprisingly low viscosity 
[1]. While the initial conditions have to be modeled [2j, the hydrodynamic evolution relies, 
at least up to first order, on the QCD equation of state (EoS) and two additional transport 
coefficients. The EoS can be obtained by first principle lattice calculations and we report on 
recents strategies of the HotQCD collaboration to calculate the EoS with the highly improved 
staggered quark (HISQ) action. Note that due to the notorious sign problem lattice QCD 
calculations are currently restricted to zero baryon chemical potential. However, the effect of a 
nonzero chemical potential can be taken into account by a controlled Taylor expansion in /U^/T 



We continue with a discussion of fiuctuations of conserved charges such as net baryon number, 
electric charge and strangeness. As one expects a critical point in the phase diagram at nonzero 
density, fluctuation observables are the most natural choice to detect critical behavior in QCD. 
Their divergent behavior can give insight in the structure of the phase diagram and has been 
proposed as a signal for the critical point [5]. Moreover, as we argue here, they can be used 
to extract the chemical freeze-out conditions as observed in heavy ion collisions [HI E]. After 
hadronization of the fireball, we expect a strongly interacting gas of hadrons with a fluctuating 
chemical decomposition. However, from a given point on the trajectory, the particle abundances 
stay flxed. This point is called chemical freeze-out. The coordinates of this point in the {T, fis)- 
diagram are usually obtained by fltting the experimentally measured particle abundances to the 
hadron resonance gas (HRG) model. Those fits have been quite successful in the past [8]. It has 
been shown that the freeze-out parameters that are obtained by HRG model fits from various 
experiments at different collision energies fall on a unique curve that can be parametrized as a 
function of the center of mass energy {y/sNN) [[Qj. 

2. The highly improved staggered quark (HISQ) action 

Lattice QCD simulations utilize a discretized version of the QCD action, which at non-zero 
values of the lattice spacing is not uniquely defined. In fact, many lattice actions exist, which 
differ in the way they handle fermion doublers, to what extent the chiral symmetry is preserved 
and also in the amount of discretization effects they possess. A widely used class of lattice actions 
is given by various improved versions of so called staggered actions. They are numerically cheap 



laii!. 



and preserve a U{1) subgroup of the SU{Nf)L x SU{Nf)ji chiral symmetry that is manifest 
in the QCD Lagrangian at zero quark masses. Results that are presented here have been 
obtained with the HISQ action [lOj, that strongly reduces lattice discretization effects over the 
naive staggered action in two perspectives. Firstly, the discretization effects stemming from 
the covariant derivative term are improved by incorporating in addition to the usual nearest 
neighbor discretization of the derivative a higher order three-link term (Naik-term). This term 
dramatically reduces discretization effects at high temperature. Secondly, the flavor (taste) 
symmetry breaking, which is inherent to staggered quark actions, is reduced by a two step 
smearing procedure, where the smeared links are projected back to the gauge group SU{3) 
after the first step. The taste breaking results in an unphysical splitting of the pion spectrum at 
given lattice spacing, which then leads to an averaged pion mass (root-mean-square mass m^^^). 
Although the lightest pion mass (Goldstone pion) is kept fixed through out our simulations at 
different lattice spacings the tti^^^ varies. It can be seen as a measure of the taste breaking 
effects and is shown in Fig. [2] (left). As one can see, among the three compared staggered actions 
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Figure 2. Root-mean-square mass of the pion spectrum as function of the lattice spacing (a) for 
the HISQ/tree, stout and asqtad action (left). Arrows indicate the lattice spacings corresponding 
to our A^T- = 12, 8 and 6 lattices (from left to right) at T = 200 MeV. The right panel shows the 
temperature ambiguity as obtained from two different scales for N^- = 12, 8 and 6 lattices. 

the HISQ action performs best. The improvement of the taste breaking is most important at 



low and intermediate temperatures and, as discused in Sec. 4.3, drastically modifies the electric 
charge fluctuations in the hadron gas phase. 

The results that are presented here have been obtained with (2+l)-flavor of HISQ fermions, 
with a strange quark mass (m^) that was tuned to its physical value and light quark masses that 
are fixed by the quark mass ratio mi/ms = 1/20. This is slightly heavier than the physical value 
mi/nis ~ 1/27. We calculate on lattices with temporal extent of A^,- = 6,8,10,12 and spacial 
extent Nu = 4Nt-. For the T = calculations needed to renormalize the EoS calculations, we 
use Nr > N„. The temperature scale was set by two independent experimentally accessible 
quantities, the Sommer scale [11] (ri) and the Kaon decay constant {Jk)- The ambiguity in the 
temperature obtained by those two quantities at finite is shown in Fig. [2] (right). 

3. Status of the EoS calculation 

The quantity most convenient to calculate on the lattice is the trace anomaly in units of the 
fourth power of the temperature Q^'^ /T^. This is given by the derivative of p/T^ with respect 
to the temperature, i.e. 



rp4 rp4 Qrp 



Since the pressure is given by the logarithm of the partition function, p/T = V^^lnZ, the 
calculation of the trace anomaly requires only the evaluation of straightforward expectation 
values. Note, however, that we have to subtract the corresponding T = result in order to 
remove the UV divergencies. I.e., for each parameter set of couplings and masses we have to 
perform an additional T = calculation, which makes the EoS determination numerically very 
demanding. All further thermodynamic quantities can be deduced from the trace anomaly, e.g. 
the pressure is obtained as 
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In Fig. [3] we show our current status of the Q'^'^{T)/T'^ calculations from A'^,- = 6,8, 10 and 12 
lattices [12\. Also shown is a continuum estimate, obtained with a different staggered type action. 
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Figure 3. Comparison of the HISQ interaction measure on = 6, 8, 10 and 12 lattices with 
the stout continuum estimate [l3] (left). The HISQ interaction measure at low (middle) and 
high (right) temperatures. Solid lines in the middle panel indicate HRG results, dashed and 
solid lines in the right panel show 1-loop and 2-loop perturbative results, respectively. 

the stout action, taken from Ref. [13]. As one can see, results obtained with the two actions 
currently are not in agreement in the temperature region 180 ^ T < 300 MeV. The discrepancies 
are most pronounced at T ~ 200 MeV, where the trace anomaly has its maximum. Note, that 
the HISQ data is not yet extrapolated to the continuum, which may reduce or eliminate the 
discrepancy. It is, however, obvious that for the pressure that is obtained by the integral in 
Eq. Q, we encounter differences of (10-15)% for T > 300 MeV between the present HISQ data 
at Nt = 12 and the stout continuum estimate [13 j . 

The middle panel in Fig. [3] shows the low temperature region of the interaction measure. 
The solid lines indicate corresponding HRG results, where from top to bottom the pion mass 
was set to its physical value and to the root-mean-square mass mj^^ corresponding to our 
Nr = 12, 10, 8 and 6 lattices. As one can see, the trace anomaly calculated in the HRG model is 
not very sensitive to the pion mass already at temperatures above 130 MeV. Furthermore, one 
sees that the lattice data approach the HRG model at T ~ 150 — 160 MeV. 

Fig. [3] (right) shows the trace anomaly obtained with the HISQ and stout actions in the high 
temperature regime. Also shown here are the perturbative results at 1-loop and 2-loop order, 
indicated by the dashed and solid lines, respectively. We find good agreement with the lattice 
data already above T > 300 MeV. 



4. Conserved charge fluctuations 

We will now continue with the discussion of fluctuations of conserved charges, which can be 
measured in experiments and calculated on the lattice. They are well suited to study critical 
behavior but can also be used to determine freeze-out conditions. 



4.I. Cumulants of conserved charge fluctuations 

On the lattice we calculate derivatives of the partition function with respect to baryon (B), 
electric charge (Q) and strangeness {S) chemical potentials, which are also known as generalized 
susceptibilities and are defined as 

{VT^) ■ xgf (T) = [^^+^+'lnZiT,^^B,^^Q,^^s)) /(a/i^a/i^Q5A|) , (3) 

with jlx = fJ^x /T and X = B,Q, S. The lattice studies are performed at fj,x = and sufficiently 
close to the thermodynamic limit. The generalized susceptibilities are intensive quantities. They 
can also be interpreted as Taylor expansion coefficients of In Z and, furthermore, be related to 
the cumulants of the fluctuations of the net charges {Nx), which are measured in heavy ion 
collisions. E.g., for the diagonal fluctuations one obtains 

{VT')-X^ = {{SNxf), (4) 

{VT^)-Xf = {idNx)'')-3{i5Nxfy , (5) 

{VT')-xi = {idNxf)-l5{{6Nx)'){i5Nxf)+30{i5Nx?)\ (6) 

with 6Nx = Nx- {Nx). In Fig. |4] (left and middle) we show our current results for the diagonal 
fluctuations of net baryon number and electric charge, obtained with the HISQ action on lattices 
with temporal extent of Nj- = 6 and 8. 
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Figure 4. Cumulants of net baryon number fluctuations (left) and net electric charge 
fluctuations (middle) up to the sixth order; quadratic cumulant of strangeness, baryon number 
and electric charge normalized by the HRG results (right). Different symbols denote different 
lattice spacings: Nr = 6 (triangles), 8 (circles) and 12 (squares). Solid lines on the right panel 
indicate spline fits to the data. For the electric charge fluctuations additional solid lines indicate 
HRG results with a modified pion mass (see text for explanation). The vertical bar indicates 
the QCD transition temperature as obtained in [T7] . 



4-2. 0(4) critical behavior 

In the chiral limit of QCD with two massless quarks the transition from the hadronic world 
to the quark gluon plasma is presumably of second order (see Fig. IT|. At vanishing chemical 



potential the critical point is expected to be in the universality class of the three dimensional 
0(4) symmetric spin model and is expected to persist even in the presence of a physical strange 
quark mass. It is illuminating to perform a universal scaling analysis that is connected with the 
critical point in the chiral limit. Here to leading order the chemical potential enters only in the 
temperature like scaling field 

* = *0 ^ ((^ - Tq)/Tq + + K^/xl + KBsfj'BlJ's) , (7) 

as a finite chemical potential does not alter chiral symmetry breaking. It is thus easy to see that 
at fix = the contribution from the singular part of the partition functions to the generalized 
susceptibilities follows the pattern [l3], 

\j.\2—n—a fo\ 
X2n ~ m , (oj 

where a is the critical exponent of the specific heat, which is small and negative. The fourth 
order cumulants will thus develop a cusp in the chiral limit, whereas the sixth order cumulants 
are divergent with amplitudes that have different signs below and above Tc. The lattice data 
is qualitatively consistent with that picture as can be seen in Fig. [4] (left and middle). A 
more detailed scaling analysis using a recent parametrization of the 0(4)-scaling function of the 
specific heat [15] is work in progress. It will allow to obtain various non- universal normalization 
constants that map QCD to the universal 0(4)-symmetric theory such as TQ,tQ,ho, kx, kxy- 
Some of them are of immediate interest, such as the transition temperature in the chiral limit Tq 
or the curvature nx that characterizes the change of the transition temperature in the direction 
of the chemical potential ^x- Recently the HotQCD Collaboration and the BNL-Bielefeld 
Collaboration performed a scaling analysis of the chiral condensate, the chiral susceptibility 
|16| [T7] as well as a mixed susceptibility . That also determined the normalization constant 

zo = h\l^^ /tQ, where /? and 6 are critical exponents, which fixes the quark mass dependence of 
the QCD transition (at ^x = 0). This was used by HotQCD in order to determine the crossover 
temperature Tc = 154(9) MeV [T7] at physical quark masses. 

4-. 3. Comparison with the Hadron Resonance Gas 

Below the QCD transition temperature we can compare our lattice results with the statistical 
Hadron Resonance Gas (HRG) model, which describes the hadronization process in heavy ion 
collisions quite successfully [8]. For the second order cumulants of baryon number, electric charge 
and strangeness fluctuations this comparison is performed in Fig. |4] (right) [19]. Here we plot 
the lattice data normalized by the corresponding HRG results. Light blue bands indicate the 
continuum extrapolations based on lattices with temporal extent Nj- = 6, 8, 12. We find that the 
continuum extrapolations of X2 ^^'^ X2 approach the HRG from below and are in agreement with 
the HRG at temperatures up to T ~ (150 — 160) MeV. The strangeness fluctuations, however, 
seem to overshoot the HRG and eventually approach the HRG at lower temperatures from above. 
Note that for the electric charge fluctuations additional solid lines at low temperature indicate 
HRG results with a modified pion mass. We have chosen the pion masses such that they agree 
with the averaged pion mass (root-mean-square mass m^^^) of the pion spectrum on the lattice 
at given A'^^. In general, we find that the modified HRG provides a good approximation to the 
electric charge fluctuations below Tc. 

A comparison of the HRG with higher order cumulants is at present not very meaningful as 
we do not have a continuum result yet. Moreover, for the electric charge fluctuations, which can 
in principle be immediately compared to the experimental results, the distorted pion spectrum 
on the lattice becomes increasingly problematic, as the higher order cumulants are increasingly 
sensitive to an increased m^^^. One thus has to use even finer lattice spacings such as Nj- = 16 
in order to control this systematic effect. In fact, we find that the HRG results for Xa with a 
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Figure 5. The leading and next-to- leading order expansion coefficients of the strangeness (left) 
and the negative of the electric charge chemical potentials (middle) versus temperature for 
r = 0.4. For si and qi the LO-bands show results for the continuum extrapolation. For S3 
and we give an estimate for continuum results (NLO bands) based on spline interpolations 
of the A'^^ = 8 data. Dashed lines at low temperature are from the HRG model and at high 
temperature from a massless, 3-flavor quark gas. The right hand panel shows NLO results for 
IJ's/fJ'B and fiQ/ fj-B as function oi for three values of the temperature. 
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4.4- Determination of freeze- out conditions 

In oder to compare the lattice QCD results with experimentally measured fluctuations one has 
to eliminate the unknown fireball volume by taking ratios of cumulants. We have recently 
proposed a detailed method to extract the remaining freeze-out parameters, i.e. the freeze-out 
temperature (T-^) and the freeze-out chemical potential (/U^g), from a comparison of lattice and 
experimental results of ratios of cumulants of the electric charge fluctuations [6^ 7J . In order to 
resemble the conditions met in heavy ion conditions we fix the strangeness and electric charge 
chemical potentials by demanding net strangeness neutrality and the correct iso-spin asymmetry: 



{Ns) = 0, {Nq) = r {Nb) . 



(9) 



Here r is a constant. We choose r = 0.4, which approximates well the situation met in Au-Au 
and Pb-Pb collisions. We fulfill the conditions given in Eq. [9] by expanding the net densities as 
well as the chemical potentials fls and fiq in terms of fiB- We then solve order by order for the 
expansion coefficients of the chemical potentials. In next-to-leading order (NLO) we obtain 



flS = Si{T)flB + S3(T)/i| + 0{fL%), 



fLQ = qi{T)fLB + qz{T)(i% + 0{fL% 



(10) 



where the results for si(T), s^{TY qi(T), 93 (T) can easily be expressed in terms of the generalized 
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(left and middle) we show those expansion coefficients. The 
upper panel shows the leading order (LO), whereas the lower panels show the ratio NLO to LO 
coefficients. The band in the upper panels indicates the continuum extrapolation based on the 
A'^,- = 6,8 and 12 data, the band in the lower panels is a continuum estimate based on A^.^ = 6 
and 8 data. We find that the NLO contributions are negligible in the high temperature region 
and below 10% in the temperature interval relevant for the analysis of freeze-out conditions, i.e., 
T ~ (160 ± 10) MeV. In fact, in this temperature range the leading order lattice QCD results 
deviate from HRG model calculations expanded to the same order by less than 15%. Note that 
one can also investigate the convergence properties of the HRG model itself. In the HRG model 
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Figure 6. The left panel shows LO (top) and NLO (bottom) expansion coefficients of i?^ 
for for r = 0.4. The bands and lines are as in Fig. gleft) .On the middle panel we plot i?^2 
versus /x^/T, including the NLO contribution, for three values of the temperature. R^^ versus 
temperature is shown on the right panel for /i^ = 0. The wider band on the data set for A',- = 8 
shows an estimate of the magnitude of NLO corrections. 



the NLO expansion reproduces the full HRG result for fig and fis to better than 1.0% for all 
values of fiB ^ 1-3. Altogether, we thus expect that the NLO truncated QCD expansion is a 
good approximation to the complete QCD results for fig and fts for /i^ < 200 MeV. 

Our results for the strangeness and electric charge chemical potentials at NLO as function of 
fj,B and T are shown in Fig. [s] (right). While fis/fJ-B varies between 0.2 and 0.3 in the interval 
150 MeV < T < 170 MeV, the absolute value oi ^q/ is an order of magnitude smaller. Both 
ratios are almost constant for hb ^ 200 MeV, which is consistent with HRG model calculations. 

We do now consider ratios of cumulants Rnm — Xn IXmi which to a large extent eliminate the 
dependence of cumulants on the freeze-out volume. After fixing and ^ls by the constraints 
given in Eq. |9j they only depend on /i^ and T. Ratios with n + m even are non-zero for ^b = 0; 
while the odd-even ratios are in leading order proportional ^b and thus vanish for fj,B = 0. 
Ratios with n + m even or odd thus provide complementary information on T-f and /u^g. The 
simplest of those ratios are given by 

Rf2 = ^ = AB(i?f/ + ^f/ Al + 0(A|j)) , (11) 
Rl ^ ^ = R^' + R^ff^l + Oif.%). (12) 

These ratios can be calculated in QCD as well as in the HRG model [14j, and eventually can 

be compared to experimental data in order to determine and //^ . Here we have evaluated 

them up to the NLO contribution 0{fl\) in the Taylor expansion of R^2 and to LO for R§^. 

For the discussion of results presented below, we concentrate on net electric charge fluctuations, 

i.e. on r!^2 and R!^i. While net baryon number fluctuations are experimentally accessible only 

through measurements of net proton number fluctuations, which may cause some difficulties 

|2H I22|. electric charge fluctuations may be easier to analyze. They are thus better suited for a 

comparison with the experimental data. 

In Fig. [g] (left) we show our results on the LO and NLO expansion coefficients of i?^. The 

bands on the lower and upper panel have the same meaning as in Fig. Isl (left). We find that 

o — 
the NLO corrections to are below 10%, which makes the LO result a good approximation 

for a large range oi jiB- Systematic errors arising from the truncation of the Taylor series for 

R^2 at NLO may again be estimated by comparing the full result in the HRG model calculation 

with the corresponding truncated results. Here we find for T = (160 it 10) MeV and jlB ^ 1-3 



that the difference is less than 1.0%. Moreover, we estimated that taste violation effects in the 
NLO calculation lead to systematic errors that are at most 5% and thus will be negligible in 
• Taylor series truncated at NLO are thus expected to give a good approximation to the 
full result for a wide range of baryon chemical potentials. 

In Fig. 6 (middle) we show the full fiB and T dependence of R^2^ including the NLO 

— o 
contribution. Obviously the ratio shows a strong sensitivity on fis but varies little with 

T in the temperature range T = (160 it 10) MeV. For the determination of (T-^,/i^) a second, 
complimentary information is needed. To this end we use the ratio R^i, which is strongly 
dependent on T but receives corrections only at 0{fi'^). The leading order result for this ratio is 
shown in Fig. [g] (right). Apparently this ratio shows a characteristic temperature dependence for 
T > 155 MeV that is quite different from that of HRG model calculations. The NLO correction 
to this ratio vanishes in the high temperature limit and at low T the HRG model also suggests 
small corrections. In fact, in the HRG model the LO contributions to R^^ differ by less than 2% 
from the exact results on the freeze-out curve for hb 200 MeV. The broader band in Fig. [6] 
(right) indicates an estimate of the NLO contribution at /is = 1 from our A^,- = 8 calculations. 

We now are in the position to extract /x^ and Tj from i?^ and R^^ which eventually will 
be measured in the beam energy scan at RHIC [23]. A large value for r'^^, i.e. R^i ^ 2 would 
suggest a low freeze-out temperature T < 155 MeV, while a value R^^ ~ 1 would suggest a 
large freeze-out temperature, T ~ 170 MeV. A value of i?^^ ~ 1.5 would correspond to T ~ 160 
MeV. A measurement of R^^^ thus suffices to determine the freeze-out temperature. In the HRG 
model parametrization of the freeze-out curve [9J the favorite value for in the beam energy 
range 200 GeV > ^/snn > 39 GeV indeed varies by less than 2 MeV and is about 165 MeV. 



At this temperature the va" 
bit, as is obvious from Fig. 
QCD at T = 165 MeV. Va: 



ues for R^^ calculated in the HRG model and in QCD differ quite a 
e] (right). While R^^ ~ 2 in the HRG model, one finds R'^^ ~ 1.2 in 
ues close to the HRG value are compatible with QCD calculations 
only for T < 157 MeV. We thus expect to either find freeze-out temperatures that are about 5% 
below HRG model results or values for R^^ that are significantly smaller than the HRG value. 
A measurement of this cumulant ratio at RHIC thus will allow to determine T-^ and probe the 
consistency with HRG model predictions. 

or any of these temperature values a comparison of an experimental value for with 
(middle) will allow to determine fi^^. To be specific, at T = 160 MeV we expect to 
= (20 - 30) MeV, if R% lies in the range 0.012 - 0.020, /x^ = (50 - 70) MeV for 

0.032 < i?^2 < 0-045 and //^ = (80 - 120) MeV for 0.05 < R^^ < 0.08. These parameter ranges 
are expected J9j to cover the regions relevant for RHIC beam energies ^/~SNN = 200 GeV, 62.4 
GeV and 39 GeV, respectively. As is evident from Fig. 6 (middle) the values for will shift to 
smaller (larger) values when turns out to be larger (smaller) than 160 MeV. A more refined 
analysis of (T'^,/i^) will become possible, once the ratios i?^ and R^-^^ have been measured 
experimentally. 
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5. Summary 

We have reviewed the current status of our equation of state calculations with the HISQ action, 
based on lattices of temporal extension A^,- = 6, 8, 10, 12. Recent results for the trace anomaly 
indicate a discrepancy of roughly 20% in the peak region (T ^ 200 MeV) between HotQCD [12j 
and the Budapest-Wuppertal results [13]. However, final conclusions can only be drawn after 
the continuum extrapolation of our data, which we will attempt as soon as the Nr = 10 and 
12 data have reached higher statistical significance. We have further discussed fiuctuations of 
conserved charges. We have shown that higher order cumulants of conserved charge fiuctuations 
become more and more sensitive to critical behavior. For the second order fiuctuations we 



have explicitly discussed their approach to the HRG in the low temperature region. Finally, we 
proposed a method to extract freeze-out parameters, i.e. the freeze-out temperature (T-^) and 
the freeze-out chemical potential (^^g), from a comparison of lattice and experimental results 
of ratios of cumulants of the electric charge fluctuations jSl [7j . At nonzero baryon density the 
method is based on a next to leading order (NLO) Taylor expansion of two different ratios. In 
general we find that the NLO contributions are below 10% in the f^s/T range relevant for BES 
at RHIC down to collision energies of ^/snn ^ 20 GeV. 



Acknowledgments 

The numerical generation of gauge field configurations has been performed on BlueGene/L 
computers at Lawrence Livermore National Laboratory (LLNL), the New York Center for 
Computational Sciences (NY-CCS) at Brookhaven National Laboratory, US Teragrid (Texas 
Advanced Computing Center), Cray XE6 at the National Energy Research Scientific Computing 
Center (NERSC), and on clusters of the USQCD collaboration in JLab and FNAL. The 
evaluations of fluctuation observables have been performed on the USQCD owned GPU-cluster 
at JLab and the GPU-cluster at the University of Bielefeld. We further acknowledge support 
by contract DE-AC02-98CH10886 with the U.S. Department of Energy, the Bundesministerium 
fiir Bildung und Forschung under grant 06BI9001, the Gesellschaft fiir Schwerionenforschung 
under grant BILAER, the Deutsche Forschungsgemeinschaft under grant GRK881 and the EU 
Integrated Infrastructure Initiative Hadron-Physics 2. 



References 

P. Romatschke and U. Romatschke, Phys. Rev. Lett. 99, 172301 (2007) ,arXiv:0706.1522] [nucl-th]]. 
[2] C. Gale, S. Jeon, B. Sclienke, P. Tribedy and R. Venugopalan, arXi v:1210.5144| [hep-ph], 
[3] C. R. AUton, et al, Phys. Rev. D 66 (2002) 074507 .hep-lat/0204010[ ! 
[4] C. Schmidt, Eur. Phys. J. C 61, 537 (2009). 

[5] M. A. Stephanov, K. Rajagopal and E. V. Shuryak, Phys. Rev. D 60, 114028 (1999) [hep-ph/9903292| . 
[6] A. Bazavov et al., arXiv: 1208. 1220 [hep-lat], to appear in Phys. Rev. Lett. 
[7] S. Mukherjee, Quark Matter 2012. 

P. Braun-Munzinger, K. Redlich, J. Stachel, In *Hwa, R.C. (ed.) et al.: Quark gluon plasma*, 491-599. 
[9] ,L Cleymans, H. Oeschler, K. Redlich and S. Wheaton, Phys. Rev. C 73, 034905 (2006). 
[10] E. FoUana et al. [HPQCD and UKQCD Collaborations], Phys. Rev. D 75, 054502 (2007). 
[11] R. Sommer, Nucl. Phys. B 411, 839 ( 1994) ;hep- lat/9310022 . 
[12] A. Bazavov [HotQCD Collaboration], 



arXiv: 1210 .6312 [hep-lat]. 

, JHEP 1011, 077 (2010), [arXiv:1007.2580l ;hep-lat]]. 



[13] S. Borsanyi et al, [BW Collaboration 

[14] S. Ejiri, F. Karsch, K. Redlich, Phys. Lett. B633, 275 (2006) 

[15] J. Engels and F. Karsch, Phys. Rev. D 85, 094506 (2012) |arXiv:1105.0584 [hep-lat 

[16] S. Ejiri et al, Phys. Rev. D 80, 094505 (2009) arXiv: 0909.5122. [hep -lat]] . 

[17] A. Bazavov et al., Phys. Rev. D 85, 054503 (2012) ar Xiv:1111.1710 ' [hep-lat]]. 

[18] O. Kaczmarek et al, Phys. Rev. D 83, 014504 (2011) arXiv:1011.3130 [hep-lat]]. 

[19] A. Bazavov, et al, Phys. Rev. D 86, 034509 (2012) arXiv: 1203.0784 [hep-lat]]. 

[20] M. M. Aggarwal et al. [STAR Collaboration], Phys. Rev. Lett. 105, 022302 (2010). 

[21] M. Kitazawa, M. Asakawa, Phy s. Rev. C 85, 021901 (2012). 

[22] A. Bzdak, V. Koch, V. Skokov, |arXiv: 1203.4529 [hep-ph]. 

[23] L. Kumar [STAR Collaboration], Quark Matter 2012. 



